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Tagged particle in single-file diffusion 


Abstract Single-file diffusion is a one-dimensional interacting infinite-particle system in which the 
order of particles never changes. An intriguing feature of single-file diffusion is that the mean-square 
displacement of a tagged particle exhibits an anomalously slow sub-diffusive growth. We study the 
full statistics of the displacement using a macroscopic fluctuation theory. For the simplest single-file 
system of impenetrable Brownian particles we compute the large deviation function and provide an 
independent verification using an exact solution based on the microscopic dynamics. For an arbitrary 
single-file system, we apply perturbation techniques and derive an explicit formula for the variance in 
terms of the transport coefficients. The same method also allows us to compute the fourth cumulant 
of the tagged particle displacement for the symmetric exclusion process. 

Keywords Single-file diffusion, Macroscopic fluctuation theory. Anomalous diffusion, Large deviations 
PACS 05.40.-a, 05.70.Ln, 05.70.Np, 05.10.Gg 


1 Introduction 


In non-equilibrium statistical mechanics, dynamical properties of interacting many-body systems play 
as important role as non-equilibrium steady states. Unlike the equilibrium properties, dynamical prop¬ 
erties depend on the initial and boundary conditions and on details of the interactions. 

The motion of an individual particle in a system of interacting particles is a fundamental dynamical 
problem in statistical mechanics even when the entire system is in equilibrium, or in a non-equilibrium 
steady state. A single Brownian particle exhibits normal diffusion—the mean-square displacement 
grows linearly with time. In a gas composed of impenetrable diffusing particles, the mobility of an 
individual particle is reduced due to the presence of other particles. In two and three dimensions, this 
effect can be significant in dense gases, but qualitatively the tagged particle still undergoes normal 
diffusion. In one dimension, however, the order of particles is preserved and transport can become 
anomalously slow. Transport in such systems is known as single-file diffusion. 

Single-file diffusion is prevalent in numerous physical, chemical and biological processes: molecular 
motion inside porous medium like zeolite , water transport inside a carbon nanotube , motion 
of tagged monomers in a polymer chain [J, sliding of proteins in a DNA sequence j^, ion channels 
through biological membranes [^, super-ionic conductors [^, etc. 

The motion of an individual (tagged) particle strongly depends on whether the dynamics is biased 
or not. A diffusive behavior emerges when the dynamics is biased [8 14 . In the unbiased case, the 
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variance of the displacement exhibits a sub-diffusive VT growth. In this paper, we focus on the 
unbiased situation^ 

Vt behavior was first derived by Harris 


17 


for the system of Brownian 


The sub-diffusive {X^) 

particles interacting only through the hard-core repulsion; the same qualitative behavior was found 
to occur 


18 19 


for the symmetric exclusion process (SEP). Since then, a large amount of work has 

Most of these studies concentrate 


30 


been devoted to the study of more general single-file systems ■ 
on the calculation of the variance, and even this requires a rather elaborate analysis when there are 
interactions in addition to the hard-core repulsion. The sub-diffusive scaling has also been verified in 
several experimental systems [3 31-35 


Additional interest has been triggered by the connection to interface fluctuations. The problem of 
tagged particle displacement can be mapped to the height fluctuation of a one-dimensional interface 
10 . The unbiased case is related to the Edwards-Wilkinson interface growth 36 , whereas the biased 


case resembles the Kardar-Parisi-Zhang interface growth 37 


An analytical treatment of single-file diffusion is a challenging many-body problem due to the 
correlations caused by the non-crossing condition. The sub-diffusive s caling of the displacement of the 
tagged particle, (A|.) ~ '/T, is easy to understand heuristically 7138 -40 . An intriguing and a bit 


counter-intuitive property of single-file diffusion is that the variance aepends on the initial state. More 
precisely, the scaling of the variance is robust, (A1|.) ^ Vt, but the pre-factor depends on whether we 
average only over the realizations of the stochastic dynamics or we additionally perform the averaging 
over initial positions [M 41 ^2 . 


Our goal is to study the higher cumulants, in fact the full statistics of the tagged particle displace¬ 
ment. The cumulants are encoded in the cumulant generating function 


A2 


A3 


Mt(A) = log = A (At), + ^ (A^+ - (A^)^ + 


3! 


( 1 ) 


as Vt, see |19[|43 
explicitly using a 


where (A^), is the cumulant and A a fugacity parameter. For example, the first cumulant is the 
average and the second cumulant is the variance. 

Fot impenetrable Brownian particles and for the SEP all the cumulants of the displacement scale 
. The same is expected to hold for more general single-file systems, as we shall show 
lydrodynamic formulation. In long time limit, the probability of the rescaled position 
of the tagged particle has a large deviation form 

= ( 2 ) 

Here is the large deviation function. The symbol x implies that the logarithms exhibit the same 
asymptotic behavior: A x B means that limT-j-oo = 1. The cumulant generating function /iT(A) 


is related to (p{x) via a Legendre transform 44, ^ 

In the present work we apply the macroscopic fluctuation theory (MET) to single-file diffusion. The 
MET was developed by Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim 46 -50 for calculating 


large deviation functions in classical diffusive systems; similar results were obtained in the context of 
shot noise in conductors 


51 52 


theoretical framework for non-equilibrium systems 


successfully applied to numerous systems 46 - 48,53 


The MET provide s a sign ificant step towards constructing a general 

Over the past decade the MET has been 
perfect agreement between the MET and 


systems 

44 

53 

46 

00 

53 

-70 
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microscopic calculations has been observed whenever results from both approaches were available. The 
MET is a powerful and versatile tool, although the analysis is involved and challenging in most cases. 

The MET allows one to probe large deviations of macroscopic quantities such as the total current. 
Intriguingly, large deviations of an individual microscopic particle in a single-file system can be captured 
by the MET 29 69 . The essential property of single-file systems—the fixed order of the particles— 

Y amenable to the MET treatment. 
Using a path integral formulation 


allows us to express the displacement of the tagged particle in a way amenable to the MET treatment. 
The MET is an outgrowth of fluctuating hydrodynamic 


71 


one associates a classical action to a particular time evolution of the system. The optimal path has 
the least action and the analysis boils down to solving the Hamilton equations corresponding to the 
least-action paths. The advantage is that within this formulation all the microscopic details of the 


If the particle dynamics is ballistic, the tagged particle exhibits a normal diffusion behavior |l5l 161. 
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system are embedded in terms of two transport coefficients: diffusivity D{p) and mobility <j{p). These 
are bulk properties of the system near equilibrium, which for an arbitrary single-file system can either 
be measured in experiment or calculated from the microscopic dynamics. 

We shall show how the cumulant generating function ^r(A), equivalently the large deviation func¬ 
tion 0(x), can be formally expressed in terms of D[p) and cr(p). The y/T scaling of the cumulants comes 
out from this formal solution. The dependence of the cumulants on the initial state is also naturally 
incorporated within this formalism. Different initial states lead to different boundary conditions, while 
the governing Hamilton equations remain the same. 

Mathematically, the Hamilton equations are a pair coupled non-linear partial differential equations 
for two scalar fields. In the general case when D[p) and <7{p) are arbitrary, the Hamilton equations are 
intractable. The only solvable case corresponds to Brownian particles where D{p) = 1 and cr(p) = 2p. 
For this single-file system we deduce a closed form formula for the cumulant generating function and the 
associated large deviation function. We analyze annealed and quenchec0 initial states. Our results show 
that the single-file system remains forever sensitive to the initial state—the large-deviation function 
has a very different behavior in the two cases. We verify the MFT predictions in the particular case of 
Brownian particles by comparing with exact results which we derive using the microscopic dynamics: 
the large deviation functions coming from these independent methods match perfectly. 

In the general case of arbitrary D[p) and (j{p) we use a series expansion method. In principle, 
the cumulants of all order can be evaluated iteratively, but in practice the calculations become very 
cumbersome as the order increases. Our analysis leads to an exact result for the variance of the tagged 
particle in terms of D{p) and <j{p). The variance in the annealed and quenched cases differ by y/2. 
This was observed for the symmetric random average process 72 and for impenetrable Brownian 


particles 42 , and it remains generally valid for single-file diffusion. Our general formula applies to an 


arbitrary single-file system and it encompasses all the results derived for specific models 17[ 18[|21 42 
as well as experimental results 33 . For the SEP, where the transport coefficients are D[p) = 1 and 
(t(p) = 2/9(1 — p), we derive an explicit formula for the fourth cumulant. Our derivation is based on the 
macroscopic MFT framework, yet we also needed exact results for the integrated current which were 
derived in Ref. [73] on the basis of an exact microscopic analysis employing a Bethe ansatz. 

The remainder of this paper is structured as follows. In Section we start from a fluctuating hy¬ 
drodynamics and show how the statistics of the tagged particle in single-file diffusion can be integrated 
into the MFT, leading to a variational problem. In Section we analyze the single-file system of 
Brownian point particles. A perturbative treatment of the Hamilton equations with arbitrary D[p) 
and a{p) is presented in Section]^ This allows us to derive a general formula for the variance of Xt- 
In Section we calculate the fourth cumulant for the SEP. We present a microscopic analysis of the 
single-file system of Brownian particles in Section Some intermediate technical steps are relegated 
to the Appendices. 


2 A hydrodynamic formulation 


The macroscopic fluctuation theory (MFT) is a deterministic re-formulation of fluctuating hydrody¬ 
namics. It has proven to be a very successful framework [50|61| for probing large deviations in diffusive 
particle systems (such as lattice gases). The MFT generalizes the Freidlin-Wentzell theory 74 of finite¬ 
dimensional dynamical systems with random perturbations to a class of stochastic infinite-dimensional 
dynamical systems. 

The simplicity of lattice gases, and generally diffusive particle systems, is that on a macroscopic level 
the full description is provided by one scalar field, the density, satisfying the diffusion equation. Taking 
into account the stochasticity of underlying microscopic dynamics leads to a more comprehensive 
description known as fluctuating hydrodynamics 71 . The basic assumption is that the behavior is still 
essentially hydrodynamic on a large (in comparison e.g. with the lattice spacing and the hopping time) 
length and time scales. In one dimension, the fluctuating hydrodynamics is based on the Langevin 
equation 

' ( 3 ) 


dtp = dx D{p)dxP + \/cr(p) ry 


^ We shall often use this concise, but a bit imprecise, description. A quenched state is a fixed initial state; 
an annealed state is actually a collection of states, say states with macroscopically uniform density, and we 
perform averaging over all these states. 
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where 77 = 'r]{x,t) is a Gaussian noise with covariance 


{ri{x, t)r](x\ t')) = S{x — x')S{t — t'). 


( 4 ) 


The diffusion coefficient D{p) and the mobility cr{p) are related to the free energy density f{p) through 


the fluctuation-dissipation relation 44 62,71 


/"(P) 


2D{P) 

o'(p) 


( 5 ) 


which is a consequence of the assumption of local equilibrium. 

The key feature of the fluctuating hydrodynamics and the MFT is that all microscopic details 
(interactions, hopping rules, etc.) are embedded into the two transport coefficients D{p) and a{p). These 
transport coefficients generically depend on the density and can be difficult to determine analytically, 
but they can be measured numerically or experimentally. 

The governing Langevin equation ([^ describes long time properties of the system where the density 
varies smoothly over the coarse graining scale. The analysis based on this equation correctly captures 
the large scale statistics of the fluctuations. For example, in the case of single-file diffusion only the 
leading time-dependence of the cumulants is correctly captured by fluctuating hydrodynamics (equiv¬ 
alently, the MFT). 

The remarkable property that allows one to apply the MFT to single-file diffusion is a simple 
relation between the position of the tagged particle and the density p{x,t) [^. Let the tagged particle 
starts at the origin at time t = 0 and moves within a time window [0, T]. Its position Xt at any time t 
is related to the density by the single-file constraint that particles do not cross each other. This gives 


>Xt[p] 


pOO 

dx p{x,t) = / dx p{x^0)j 

Jo 


( 6 ) 


thereby defining the tagged particle positioi|^as a functional Xt[p] of the density profile p{x, t). Equiv¬ 
alently, we can write 

rXt[p\ roo 

/ dx p{x,t) = / dx (^p{x,t) — p{x,0)) . (7) 

Jo Jo 

In contrast to (§, both integrals in Q are convergent. 


2.1 Variational formulation 

The generating function of the cumulants of the tagged particle position Xt at time T can be written 
as a path integral of the density profile as 

= (^jV[p\ 5{dtp-d,[D{p)d,p+,/^) 77 ]( 8 ) 

The Dirac delta function enforces the validity of the Langevin equation ([^ . The angular brackets denote 
the averaging that can consist of two parts: over the initial density field p{x, 0 ) and over realizations of 
the stochastic noise ri{x,t) in the time window [0,T] (i.e., the averaging over the history of evolution 
of the density profile). 

If we allow the initial state to fluctuate, we must include the probability Prob[p(a:, 0)] of the initial 
density profile. It is useful to define the function 

F[p(x, 0)] = - log (Prob[p(a;, 0)]). (9) 

For an initial state at equilibrium (annealed case), F[p\ is related to the free energy. For a quenched 
initial state, there are no fluctuations of the initial profile at t = 0 and we take the initial density 
profile to be the uniform profile with density p. 

^ Xt[p] is not uniquely fixed when there are regions with density equal to zero, but such configurations are 
highly improbable. 
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The following analysis is essentially the Martin-Siggia-Rose formalism 75 76 . Incorporating the 
average over the initial state, the generating function becomes 

= J V[p] (^s(^dtp - d,[D{p)d,p + ./^) p] 

The subscript p denotes the average over the history of noise within the window [0, T] and F appears 
only in the annealed case. Replacing the delta function by an integral over a field p we get 


= / 2?[P] e 


AXT[p]-F[p(a;,0)] 


J VIp] e~ -^0 '^*-^-°^dxp{9tp-d^[D(p}d^p+^/a{p) 17]) 


The average with respect to the Gaussian variable p{x, t) is computed to yield 

= [ V[p,p] 


with action 


St[p,p] 


and Hamiltonian 


-XXt[p]+F[p{x,0)] + 
H{p,p) = '^^{dxpf 



pdtp-H{p,p) 


-D{p) [dxp) {dxp). 


( 10 ) 

( 11 ) 


(In deriving these equations, we assumed that p{x,t) —>■ 0 when x ±oo.) 

Let us rescale time by T and space by VT. The tagged particle position is also rescaled by a/T, 
so the action St[p, p] is proportional to a/T and we can write St[p,p] = VTS[p,q]. This simple 
observation leads to the anomalous scaling in single-file diffusion. Indeed, for large T, the path integral 
is dominated by the least action. The cumulant generating function is then given by 


Pt{X) = -VTS[p,q ], 


( 12 ) 


where we denote by (p, q) = (p, p) the optimal paths of the least action]^ Equation (12) implies that 
at large time all the cumulants of the tagged particle position scale as \/T. Note that the analysis 
captures only the leading T dependence of the cumulants, and the sub-leading terms of the cumulants 
come from the correction to the saddle point approximation. 


Paths of least action 


To determine the paths of least action we consider a small variation around {p,q) as p = q + Sp and 
p = p + 6p. The variation of the action is then 

pT noo / r'^ /*oo / 

SSt = J dt J dx\dtq — —\5p{x,t) — j dt J dx (dtp + \ S p{x, t) 


da; I — A 
\ 

dxl — X 


SXt SF 
dq{x,0) 6q{x,0) 

dXj' 


-p(a:,0) (5p(a:,0) 


Sq{x,T) 


+ p{x,T) ]6p{x,T). 


(13) 


For the action to be stationary, 6ST[p,q] = 0, integrals in (13) must vanish. Since Sp{x,t) and 5p{x,t) 
are arbitrary, vanishing of the integrals in the first line in (13) lead to the Hamilton equations 


Otq = -j- 

op 


and 


dtp = 

6q 


(14) 


4 


The choice of notation, p and q, hints at the Hamiltonian nature of the governing equations (141. 
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The analysis of the boundary terms gives the boundary conditions for the optimal fields {p, q). The 
density Sp{x, T) at final time T is unconstrained, so the vanishing of the integral in the third line in 
(13) leads to 

SXt 


p{x,T) = X 


Sq{x,T)' 


(15) 


The vanishing of the integral in the second line in (13) depends on the way the system is prepared. In 


the quenched case, 6p{x, 0) = 0 by definition so the integral vanishes. In the annealed setting, Sp{x, 0) 
is arbitrary, therefore its pre-factor must vanish. Thus 


g(a;,0) = p 

, , , SXt SF 


quenched case, 
annealed case. 


(16) 

(17) 


The first equation states that the initial density in the quenched case is assumed to be uniform. 

For the rest of the paper, we denote the position of the tagged particle corresponding to the least 
action path by 

Y = XT[q]. 

Note that Xt is in general a random variable depending on the history of the tagged particle, but Y 
is a deterministic quantity. 

Overall, the stochastic problem of characterizing statistics of the tagged particle position reduces 
to solving a variation problem. We now write the Hamilton equations and the corresponding bound¬ 
ary conditions by computing the functional derivatives. Different boundary conditions emerge for the 
quenched and annealed settings. 


Quenched case 

The Hamilton equations read 


dtp +D{q)da;xP =-^-^{dxp)^ (18) 

dtq - dx{D{q)dxq) = -dx{cr{q)dxp) (19) 

where (T'{q) = da{q)/dq. Using we compute the functional derivatives of Y (see Appendix 

and obtain the boundary conditions 

q{x,0) = p and p{x,T) =e{x - Y), (20) 


where &{x) is the Heaviside step function. 

The corresponding minimal action (10) yields the cumulant generating function 


Me(^) ^ ~ J j 


dx 


pdtq - ^{9xp) + D{q){dxq){dxp) 


Hereinafter we use the subscript Q to denote the quenched case. (For the sake of clarity, the subscript T 
denoting the time variable is omitted from the following formulas.) The cumulant generating function 
can be further simplified by using Eqs. (18)-(19) and integrating by parts. One gets 


PQiX) = XY- 


dt 


dx'^^(dxp)^ ■ 


( 21 ) 
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Annealed case 


The Hamilton equations are the same as in the quenched case, Eqs. (18)-(191. The boundary conditions 
are different. To derive them we use the well-known relation (see 44|49 and showing the consistency 
with the MFT and the fluctuation dissipation relation) 


F[p(a;,0)] = 


pp{x,o) on(r] 

dx / dr —(p(a;, 0) — r). 

. Jp ct'-' 


(r) 


Computing the functional derivative of T and 0)], we recast (15) and ( [17| ) into 

A 


pix,T) = 
pix,0) = 


qiy,T) 

A 

QiY,T) 


e>(x-Y), 
e>(x) - 


i(x,o) 2D(r) 
dr—— . 
a{r) 


( 22 ) 

(23) 

(24) 


function for the annealed initial state becomes 
PaW = XY- 


The least action (10) can be simplified using (22) and the Hamilton equations. The cumulant generating 


r4.r\r^-PI4(,(..0)-r)- r,tr 

-OO Jp cr(r) V J Jo J-OC 


dx- 


(d.p) ■ (25) 


Hereinafter the subscript A denotes the annealed case. 


A Remark on Symmetry. It is instructive to look at the symmetry properties of the action and optimal 
fields. One can notice symmetry relations 


q-x{-x,t) = q\{x,t) 
p-x{-x,t) =p\{x,t) - 

Y-x = -Yx 
Pt{—X) = mt(A) 


A 


qx{Yx,T) 


(26) 

(27) 

(28) 
(29) 


underlying the variational problem in both quenched and annealed settings. Since /i(A) is an even 
function, all odd cumulants of the tagged particle position are zero. This is the consequence of the fact 
that the microscopic dynamics is unbiased. 


3 Brownian particles with hard-core repulsion 

It is hard to analyze the variational problem for a general single-file system characterized by density- 
dependent transport coefficients D{p) and (t(p). The only system that is amenable to analyses, both 
macroscopic (based on the MFT) and microscopic, is the system of Brownian point particles with 
hard-core repulsion (Fig. [^. The system was first studied by Harris 17 who used a mapping to non¬ 
interacting particles to compute the variance of the tagged particle position. The problem was later 
studied by many other authors (see e.g. ^25 42,77 ). 

For Brownian particles the diffusion coefficient is constant, D{p) = 1, and the mobility is a linear 
function of density, a{p) = 2p. The Hamilton equations 19) become 


dtP + dxxP= -{dxpY, 

^tq Sxxq — Ux(j2qUxP^ • 


(30) 

(31) 


A canonical Hopf-Cole transformation, {p,q) —>■ {P,Q) = {e^,qe~^), is known to simplify 

Eqs. ([30|)-([M]). The new conjugate variables satisfy the Hamilton equations with H = —{dxP){dxQ). 
In the (P, Q) variables, the governing equations are (linear) anti-diffusion and diffusion equations: 


dtP + dxxP = 0 and dtQ - dxxQ = 0. 


(32) 






























Fig. 1 A sample trajectory of Brownian point particles with hard-core repulsion. The trajectories may come 
infinitely close but never cross each other, keeping the order of particles unchanged. 


Solving these equations and returning to the original variables {p, q) we arrive at a formal solution 


p{x,t) = log 


dz e 




exp 


4(T-t )) 


y/47r(r-t) 


/ oo 

dz g(z,0) 

-OO 


/ (z— 


y/Airt 


(33) 

(34) 


applicable to both quenched and annealed settings and to arbitrary p(x,T) and q{x,0). 


3.1 Quenched case 


Let us re-write (20) as 

q{x,0) = p and p{x,T) = B 0{x — Y). 
Here we have used a shorthand notation 


qiy,T)- 

The quantity Y is related to q{x,t) via 0 which becomes 

'q{x,T) 


pOO 

y ~ j 


- 1 


(35) 

(36) 


(37) 


in the quenched case. 

The boundary conditions ( [3^ depend on the solution itself. Let us proceed by treating B and Y 
as parameters to be determined later. The optimal fields in terms of B are 


p{x,t) = log 


1 + (e^ - 1) Jerfc ( J ^ 



(38) 


(39) 
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where erfc(a;) is the complementary error function. Substituting D{p) = 1 and <j{p) = 2p into the 
formula for the cumulant generating function (21) leads to 

pT poo 

Pq{X) = XY — dt dx q{dxpY ■ 


This expression for pq{X) can be simplified thanks to identity 

qidxpf' = dt{q p) - dx {pdxq - qdxP - 2qpdxp ), 
which results from Eqs. (30|-(31). Using (40) and p{x,t) —>• 0 as a; —>■ ±oo, we obtain 

/ OO pOC) 

dx q{x,T)p{x,T) + / dx q{x,0)p{x,0). 

-OO j —OO 


Plugging into this formula p(a;, T) and q(x,0) from (35) we get 

pOC) pOO 

Pq{X) = X Y — B / dxq{x,T) + p / dx p{x,0). 


(40) 


It is clear from Eqs. (38)-(39) that the functions q{x, T) and p{x, 0) approach to non-zero constants as 
X ^ OO. This implies that the integrals in the above formula are not convergent although their linear 
combination is well defined. To write P-q{X) in terms of convergent integrals we use the definition of 
Y in (37) and subsequently rearrange the integrals and obtain 

pOO pOO 

Pq{X) = (A — pB) Y + p J dx {p{x, 0) — B) + p J dx p{—x, 0). 

The integrals are now written in terms of p{x, 0) and, using (38), the formula for the cumulant gener¬ 
ating function becomes 


Mc(A) 

Vir 


= {X-pB)y + p d^og 1 -k 


o-B 


- 1 


erfc(5) 


+ p d^og 1 + 


e" - 1 


erfc(0 ■ 


where y = Y and ^ = xjsf^. This is further simplified by using erfc(^) -I- erfc(—= 2 to give 


dQ(A) 

V4T 


= (A — pB) y + p I d^ log ( 1 -I- sinh^ ( ^ ) erfc(^) erfc(—) . 


(41) 


The reason for writing the formula in this form will become clear shortly. 

So far, Y and B have been treated as parameters. One relation between these parameters is obtained 
by inserting ([39|) into (ISTl): 


y = 


e ^-1 


de 


erfc(^ - y) erfc(?/ - ^ 


2 -I- (e^ - 1) erfc{y - ^)' 

Changing the variable ^ t C + y and using erfc(^) -|-erfc(—^) = 2 we transform the above formula into 

e~ 


y = 


— I erfc(^) erfc(—^) e ^ —1 erfc(^) erfc(—^) 


dC 


/o 


2 -I- (e-® — 1) erfc(^) 


d^ 


/o 


2 -I- (e“^ — 1) erfc(^) 


Massaging this formula one arrives at a more neat form 


d = ^ ^ d^ log (^1 + sinh^ erfc(0 erfc(-0^ ■ 


(42) 


A similar self-consistent way of determining B using relation ( |36[ ) does not lead to a unique value for 
B. This is because the solution q{x,t) in (39) is singular at {x,t) = {Y,T). A graphical representation 
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Fig. 2 The optimu m d ensity field q{x,t), Eq. (39 1 , at different times. We set p = 1, B = 2 and T = 1, and 
determined Y from \42\ . The density at initial time t = 0 starts from the quenched uniform profile q{x,0) — 1 
and as time approaches t = T, the profile develops a sharp jump at the position of the tagged particle Y ~ 1.54, 
indicating a discontinuity of the function q{x,T) at x = Y. 


of this singularity is shown in Figure The problem is analogous to the case of a diffusion equation 
with a step initial profile: The solution at any time t > 0 is independent of the precise value of the 
initial profile at the position of the step. 

For any value of B, Eqs. (38)-(39) give a solution of the Hamilton equations. Only one solution 
corresponds to the minimum action. This solution can be determined by optimizing the action with 
respect to B, i.e., by imposing 


dB 


= 0 . 


Combining this with (41) and (42) we find the optimal B: 

A 

P ' 


B = 


Putting above in (41) we arrive at an explicit formula for the cumulant generating function: 

pqW 




= P 


log 1 + sinh 


^ ) erfc(0 erfc(-0 


(44) 


(45) 


The large deviation function 


The large deviation function is related to Pq{X) via the Legendre transform: 4’Q{y) = sup_), y\y 
Using (45) and (42) we can represent the large deviation function in the parametric form 

Ve(y) (^) erfc(-0^ 

with B determined from the optimality requirement 

d(j)Q 


dB 


= 0 . 


%/4T ) 


(46) 


(47) 


An equivalent representation of (f>Q(y) is 




df log 


'-V 




erfc(^) 


d^ log 




(48) 


All previous results have been derived using a macroscopic approach. In Section we show that 
the same expression for (f>Q follows from an exact microscopic analysis. This is reassuring since the 
macroscopic approach is not fully rigorous, yet much more widely applicable than exact analyses which 
are limited to simplest systems. 
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3.2 Annealed case 


The major difference with the quenched case comes from the boundary conditions (24|-([23|). When 
D[p) = 1 and a{p) = 2p, the boundary conditions become 


q{x, 0) = p exp {p{x, 0) — B 0{x)), 
pix,T) = B e{x-Y). 


The parameter B is again defined in (36) and Y is determined from Q , equivalently 

pY poo 

/ dx q{x,T) = / dx {q{x^T) — q{x,Q)'). 

Jo Jo 


(49) 

(50) 


(51) 


The boundary condition depends on the solution itself and similar to the quenched case, a solution of 
the optimal fields is found by treating Y and B as parameters. They are later determined using the 
solution. Substituting the boundary condition in the general solution (33|-(34) leads to a formula of 
the optimal fields in terms of the Y and B. 


p{x, t) = log 1 + (e® - l) ^erfc j 


F-; 


y4(T-t) 

X — Y 






(52) 

(53) 


Equation (25) for the cumulant generating function becomes 

PaW = X Y — f dx f 
J —oo J o 


rq(oo,0) 


dr 


q{x,0) 


-1 - 


f°° 

dt / dx qipxpY 

J —OO 


where we have used again D{p) = 1 and cr(p) = 2p. This formula for paW can be rewritten, thanks 
to the identity (40), as 


Pa{X) = XY— / dx g(a;, 0) log 


' —OO 

poo 


q{x,0) 


dx{q{x^ 0) — p) 


/ OO 

dx q{x, 0)p{x, 0), 

-OO 


(we have also taken into account that p{x,t) = 0 at x —t ±oo at all time t). Using (49)-(51| one can 
greatly simplify the above expression; 


/ OO 

dx (q(x, 0) - p). 

-OO 


Combining this with (|53|) we find 


= Xy + p 


Pa{X) 

Vtt 


- 1 


d^ erfc(^) + p 


e -^-1 


d^ erfc(^). 


So far, Y and B were treated as parameters. Using (53) one can recast (51) into 

POO pOO 

2y = (e^ “ l) / erfc(0 - {e~^ ~^) erfc(0 . 

J y J — y 


(54) 


(55) 


The second parameter B cannot be obtained by evaluating g(x, T) at x = F because this function 
is singular. (This singularity is evident from Figure]^) As in the quenched case, the parameter B has 
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Fig. 3 The optimum density field ci(x. t) in (531 for the annealed case at different times. We set Y = 1.5 and 
T = 1 and determined B from Eq. ( |55[ ). As time approaches t = T = 1, the profile develops a sharp jump at 
the position Y of the tagged particle, indicating a discontinuity of the fun ctio n q(x,T). Note that the profile 
q{x, t) is symmetric under {x, t) ^ {Y — x,T — t) as can be seen from Eq. (|53[). 


to be determined from an optimization criterion = 0. Using this together with (54) we derive a 
relation between B, Y and the fugacity parameter A: 


p \ 

— = (e^ — l) erfc(y) — (e”'® — l) erfc(—?/). 


To derive (56) we used 


erfc(^) = e 


-B 


erfc(^) 


(56) 


(57) 


which follows from (55) in conjunction with the elementary identity 

poo poo 

/ d^ erfc(0 - d^ erfc(0 = 2y. 


(58) 


Equations (54)-(56) constitute a parametric solution for the cumulant generating function in the 
annealed setup. Unlike the quenched case, we don’t have an explicit formula for /r^(A). Another 
parametric representation 


dAW 


A 1 — e 


B 


p 1 + e 


B 


= 1 + 2y 


P 


L-Jy 

„-B\ 


y, 

du erfc(u) 


erfc(y) 


has been reported in Refs. 29 77 . It can be obtained from (54)-(56l using (57) and (58). 


The large deviation function 




The large deviation function is again the Legendre transform: (fAiv) = sup;, ^Ay- 

it with (54) we obtain 


^. Combining 


2(('^(y) 


pOO pOO 

Jy J-y 
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Using Eq. (571 we eliminate the dependence on B and arrive at the following explicit formula for the 


large deviation function 




erfc(5) — 


erfc(^) 


(59) 


This result will be re-derived in Section using an exact microscopic analysis. 


3.3 Comparing annealed and quenched settings 
The cumulants 

By definition Q, the cumulants are obtained by expanding the cumulant generating function. For the 
quenched initial condition this expansion of AiQ(A) is simple to generate from the explicit formula (45). 
We write below the first few cumulants of the tagged particle position, 


{xl)^ = -hVf, 

\ I c p 

= ^ h- h) 


Vt, 


(60) 

(61) 

(62) 


where dz [erfc(z) erfc(—z)]*^. The first two integrals Ii and I 2 are known 79 leading to 




- 1 




(63) 

(64) 


For the annealed case, the generating function has a parametric form. To make a series expansion 
in powers of A, we first use Eq. (55) to eliminate B from (54) and subsequently make an expansion of 


the /ryi(A) in powers of y. Then we use (55) and (56) to make an expansion of y in terms of A leading 
to an expansion of (J-aW in powers of A. The first three non-trivial cumulants are: 




' Vt, 


pV^ 




= y ( 


408 180 


-bl8 


--Vt. 


(65) 

( 66 ) 

(67) 


The large deviation function 

The large deviation function depends on the setting. Strikingly different asymptotic behaviors emerge 
in the annealed and quenched settings (Figs. and [^. In both cases, the large deviation function has 
a non-Gaussian tail. More precisely, the large deviation function (f>{y) grows linearly in the annealed 
case and cubically in the quenched case (see Fig. [^: 


My) - p \y\, 


My) ^ I lyl^- 


( 68 ) 
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- 2-10 1 2 


y 

Fig. 4 The large deviation function (j>{y) of the tagged particle position in the annealed and qnenched settings. 
In both cases, the initial average density is chosen to be uniform p = 1. Note that the quenched large deviation 
function is always larger, see Appendix for a theoretical explanation. 




Log[y] 

Fig. 5 Asymptotics of the large deviation functions for large values of y. The dashed straight lines correspond 
to the power laws y^ and y, for the quenched and the annealed case, respectively. 


In the annealed case, the large y is derived from (591 and the asymptotic relations: 

1 


dz-eric{z) = e 


0{y-^) 


4\/^ 


In the quenched case, the large deviation function has a parametric form, making the analysis a bit 
more involved. We use (461 and (47) and notice that large y corresponds to large B. The integral in 

Vb 


(46) grows as for large B. Indeed, we divide the range of integration, dx = dx + J^dx, 

use the X —>■ oo asymptotic 


erfc(x) erfc(—x) = e 


I— +0 ^ ^ 


y/TT X 


J dx log ^1 + sinh^ erfc(x) erfc(—x)^ = J dx [B — x^) + 0{Vb) 


and write 
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and find This is the dominant contribution. The second integral contributes terms of order \ jB 

only. Thus p~^4>Q{y) — By — which in conjunction with (471 leads to the asymptotic in ([6^. 


A similar asymptotic dependence is also observed in the large deviation function for the current in 
the symmetric exclusion process with step initial condition 


59 It would be interesting to study 


how universal are these power-law tails. An analogy between the displacement of the tagged particle 


and current suggests that the tails are non-universal, as for the extreme current fluctuations 63 64 


4 Calculation of the variance for general single-file systems 


For a single-file system with arbitrary I?(p) and (t(p) it is impossible to solve the Hamilton equations. 
One can try to seek an asymptotic solution [62] as a series expansion in powers of A. This series solution 
is then used to expand the cumulant generating function. 

The expansions for the optimal fields read 


q[x,t) = p + \ qiix,t) + q2{x,t) + 

p{x,t) = A pi{x,t) + A^ P2ix,t) H-, 


(69) 

(70) 


In the zeroth order, the solution is deterministic: q{x,t) = p and p{x,t) = 0. This follows from the 
governing Eqs. ([T8|)-([T^ both in the quenched and annealed cases. In the first order 


dtPi + D{p)d^^pi = 0, 

dtqi - D{p)dj;j;qi =-(T{p)dxxPi- 


Using (21) and (25) we see that the variance 


dt / dx{dxPif 


(71) 

(72) 


(73) 


explicitly depends on pi,Fi,F 2 ; it also on qi (through Yi,F 2 ). Since Y is an odd function of A, see 
(28), the zeroth order term vanishes while the first order term is given by 


y, = i 


dx {qi{x,T) -qi{x,0)) 


(74) 


as it follows from (6j7). It remains to determine Fjj- In the quenched case, F does not appear (one can 
formally set F = 0 in the quenched setting). In the annealed case we use (22) and find 


F, = 


Djp) r 

<^{p) J-c 


dx{qi{x,0)) . 


(75) 


We now determine the variance (X^) by solving (0-(l^ with the specific boundary conditions 
for the quenched and the annealed initial states. 


4.1 The quenched initial state 


Plugging the series expansion into the boundary condition (20) yields 

Pi{x,T) = p~^0{x) and (?i(x, 0) = 0 


(76) 


in the linear order. With the above boundary condition the solution of 0 can be expressed as 

dxPi{x, t) = p-^ g(0, T\x, t), (77) 

where g is the diffusion propagator 


g{z,T\x,t) = 


y4^F(p)(F-t) L 4F(p)(F-t)_ 


{z-xf 


(78) 
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for 0 < i < T. Since Pi{x, f) = 0 at a: —)■ —oo, we obtain 


Pi{x, t) = - [ dz g{0, T\z, t) = ^ erfc ( ■ ' | . 

pj-oo ^ 2p \2^Dip) (T-t)] 


(79) 


Taking into account the solution for qi{x,t) in ( [7^ can be written as qi{x,t) = —d^ipix^t) with 

(80) 


tlj{x,t) = f dr [ dz g{x,t\z,T) g{0,T\z,T). 


P Jo 


These are the only two quantities required for simplifying (X^) given by Eq. (731. Recalling that 


( 7 i(a;, 0 ) = 0 (and that formally F 2 = 0 in the quenched initial state) we obtain 

o'(p) 


^(^T)e = ^y dx qi{x,T) - J dt J dx{d^pif . 


Using the solution for qi{x,t), the first integral yields 


1 

- dx qi{x,T) = 
P Jo 


^T) 

P 


(81) 


(82) 


where we have taken into account that iJix, T) vanishes at a; —)■ 00 . The second integral in (811 becomes 

r f°° V’( 0 ,T) 


f 2 

/ dt dx[dxPi) = 
Jo J — 00 


2p 


(83) 


Combining these two results we reduce the expression for the variance to 


1 




2 ! 


^ (j{p) 

2p 2p2 


’ nOO 

dt / dx [{g{0,T\x,t)Y 
J —00 


We use (78) and compute the integral. The final expression for the variance, 

/y2\ _ <^{P) Vr 

^ P^ yra(p)’ 

is valid in the quenched setting in the general case of arbitrary D[p) and cr(p). 


(84) 


4.2 The annealed initial state 


Inserting the series expansion (69)-(70) into the boundary condition (24)-(23l we obtain the boundary 
conditions in the first order 


Pi{x,T) = 


and gi(a:,0) = {pi{x,0) - Pi{.x,T)]). 

p 2D{p) 


The variance (731 becomes 

[ {qi{x,T) - qi{x,0))dx - [ dx{qi{x,0)y 

^ PJo cr(p) 

/ \ pT poo 
a{p) 


dt / dx{d^pi{x,t)) , 


where we have used (741 and (75). 


(85) 


( 86 ) 


The following the analysis is interwoven with the quenched case in a way that allows us to avoid 
the explicit computation of most integrals. First we note that in both cases the equation for pi ( x, t) 
and the boundary condition on pi{x,T) are identical. Therefore (79) remains valid. The equation (72) 
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for qi(x,t) is also the same in both cases, but the boundary condition (85) is different. The governing 

(87) 


equation ( [72[ ) is linear in qi{x,t), so we can write the solution as a sum 

qi {x,t) = qi {x,t) + qh {x,t), 

where qi{x,t) is the solution of the inhomogeneous equation 

dtqi - D{p)dxxqi = -cr{p)dxxPi with q^(x, 0) = 0, 
and qh{x,t) is the solution of the homogeneous equation 

dtqh - D{p)dxxqh = 0 with qhix,0) = ^^^\ (pi(x,0) - pi{x,T)). 

2D[p) 


( 88 ) 


(89) 


Comparing with the quenched case we notice that qi(x,t) is same as qi{x,t) in the quenched case: 
qi{x,t) = —dxip{x,t) with ’tp{x,t) given by (80). Further, using (81) and ( [86| we find that the variance 
for the annealed case is related to the variance in the quenched case via 


1 /^2\ ^ /w2\ _ 


^(p) J-c 


1 

dx{qhix,0)y + - dx {qhix,T) - qh{x,0)). (90) 

. P Jo 


The difference depends only on the homogeneous solution qh(x,t). 

The second term on the right-hand side of Eq. (90) can be simplified using an identity 


iy dx[qh{x,T) - qh{x,0)) = J dx{qh{x,0)Y . 


(91) 


This identity is proved by using (89), then noting that dxqh{x,t) pi{x,t) does not depend on 
time because pi and qh satisfy adjoint equations, and finally that recalling pi{x,T) = 0(x)/p\ 


mp) r 

C^(p) J-c 


dx[qh{x,0)Y = 


' —OO 

poo 


dx qh{x,0){pi{x,0) -pi{x,T)) 
dxpi{x,T){qh{x,T) - qh{x,0)). 


Thus 




D{p) r 

^(P) J-c 

^{p) r 

mp) J-o 


dx{qh{x,0)Y 

) 

O 

dx{pi{x,T) -pi{x,0)y 


(92) 


where in the last stm we have used the boundary condition from Eq. (89). The last integral is computed 
using (79) to yielcjj 

a{p) f V2-l' 


1 /^2\ 1 /^2\ ^ip) (' 




D{py 


Combining this with (84) we establish a very simple general relation 


"t)q 


(93) 


(94) 


between the two initial states. This relation and Eq. (|84| are applicable for any single-file diffusion, i.e., 
for arbitrary transport coefficients|^To verify the range of applicability of Eq. (84) let us compare with 
earlier work. In the simplest case of Brownian particles we recover the well-known expression for the 


variance (for more about Brownian particles see Refs. 17,20 25 ^). For the SEP we also reproduce 


® The difference — (X|.)g is positive, as shown in Appendix 

® Our derivations are based on the MFT which assnmes local equilibrium. This is expected to be correct for 
systems with short-range inter-particle interactions. 
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the well-known result 18 . Let us also compare with the variance in single-file colloidal systems studied 

(95) 


in Ref. 21 . For comparison, we need the relation 

a{p)=2S{Q,Q)D{p) 


where S'(0, 0) is the structure factor [^. With this the variance for the annealed case becomes 

25(0, 


= 


A p2 


(96) 


which is identical to the one presented in 21 . For experimental measurements, it is convenient to 
express the variance in terms of the isothermal compressibility k. Since 5(0, 0) = I3~^k p^, 


{xl)^ = 2 /)-■ » VT. 


(97) 


where /3 is the inverse temperature. 


5 Fourth cumulant of a tagged particle in the SEP 

The symmetric exclusion process (SEP) is a diffusive lattice gas. Each site is occupied by at most one 
particle and each particle attempts to hop to neighboring empty sites with unit rate. For the SEP, the 
diffusion coefficient is constant, D{p) = 1, while the mobility has the symmetric form cr(p) = 2p(l — p) 
reflecting the p 1 — p mirror symmetry of the SEP. The formula for (t{p) can be derived using the 
fluctuation-dissipation relation (see Appendix [C|). 

We employ a perturbative approach and compute the fourth cumulant of a tagged particle position. 
It proves useful to consider a one-parameter class of models with D{p) = 1 and a{p) = 2p(l — ap). 
The parameter a varies in the range 0 < a < 1, so the mobility remains positive. The case of a = 0 
describes the single-file system of Brownian particles, whereas a = 1 corresponds to the SEP. 

Mathematically, we want to solve equations (18)-(19). The corresponding boundary conditions 
depend on the solution itself. As discussed in Section a way to solve the equations is by treating 
B = X/q(Y,T) as a parameter. This determines the optimal fields in terms of B which we denote 
by qB{x,t) and pBix,t). The solution does not explicitly involve the fugacity parameter A which 
only appears through the dependence of i? on A. A straightforward implementation of the definition 
of B does not lead to a unique value, as the function q{x,T) is singular at {x,t) = {Y,T). The value 
corresponding to the least action path is obtained by an optimization condition ^ = 0. This expresses 
B in terms of A and leads to a parametric formula of the cumulant generating function 


p = XY{B) - R{B) 


with 


A 


dY{B) 

dB 


dR{B) 

dB 


(98) 


The function R{B) is defined in terms of the optimal fields qB{x,t) and pB{x,t), Eq. (21) in the 
quenched case and Eq. ( |25[ ) in the annealed case. A detailed implementation of this procedure has 
been presented in Section for the system of Brownian particles where an exact solution was possible 
as the corresponding a(p) is linear. 

To proceed with the analysis for the quadratic cr(p), it is instructive to recall the symmetry prop¬ 
erties of the optimal fields ([M)-(p^: 


q-B{-x,t) = qB{x,t) and p-B{-x,t) = pB{x,t) - B, 

whereas Y{B) and X{B) are odd functions of B. Combining all together shows that R{B) is an even 
function of B. This is consistent with the fact that the dynamics is unbiased. 

It proves convenient to use B as & primary expansion parameter; at the end one can re-expand the 
results in terms of A and compute the cumulants. The aforementioned symmetry properties allow us 
to seek Y{B), R{B) and X{B) as the following expansions 

YiB) = YiB + YsB^ + 0{B^) 

R{B) = R2B^ + RiB'^ + 0{B^) 

A(B) = AiB-f AaB^-t 0(5^). 


(99) 
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The last formula is also equivalent to 


Ai 




( 100 ) 


Substituting the series expansions in the first equation of (98) leads to a formula for the cumulant 
generating function in powers of B, 

= (AiFi - + (A1F3 + A3T1 -Ra)B^ + -- - . 

-1 


dR fdY 


The second equation of (98), rewritten as A = -— 

— dB V dB 


Ai — 


2i?2 


and 


A3 — 


leads to 
4i?4 - 3Aiy3 


Yi Yi 

To determine the cumulants we use the above in the series expansion 0 and obtain 

(i W), = )i(-/i4 + Air3)^ 

These formal expressions hold in both annealed and quenched cases. 

The boundary condition for the optimal fields derived in (24) and 

rqsixfi) 2D(r) 

Pb(x, 0) = B 0{x) + dr —= B 0(x) + In 

Jp <^ir) 

Pb{x,T) = B 0{x-Y), 


( 101 ) 

( 102 ) 

(103) 


can be rewritten as 

qB{x,0){l - ap) 
p[l - aqB{x,G)) 


(104) 

(105) 


where we used the definition of B in (36). To proceed, we write an expansion of the optimal fields in 
powers of B as 


qB = P + qiB + 92.6 + qsB'^ + 
Pb =PiB + P2B‘^ + psB^ H-. 


(106) 

(107) 


Note that the qk{x,t) and pk{x,t) in the above formulas are different from those in (69)-(70) which 
were obtained using expansions in powers of A, and not B. 

A straightforward computation of the solution to different orders is tedious, it involves difficult 
integrals. We circumvent this by drawing comparison with the a = 0 case where an exact solution is 
available. We first illustrate this trick by computing the solution in the linear order in B. This will 
give us the second cumulant which we can compare with already known results (which were derived in 
the previous section in the general setup). Then we shall compute the forth cumulant. 


The second cumulant 


The governing equations in the first order in B are 

{dt + d^^)pi = 0 , 

{dt - dxx)qi = -2p(l - ap)dxxPi- 

To determine the corresponding boundary conditions we use a formal expansion of the step function, 

0{x -Y) = 0{x) - 6{x)Y + ^6'{x)Y^ + ■■■ . (108) 

The boundary conditions are 

Pi{x,T) = 0{x) and qi{x,0) = p [l - ap) {pi{x,0) - pi{x,T)). (109) 
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The solutions to the governing equations are almost the same as in the Brownian case: 

pi{x,t) =pi(x,t) and qi{x,t) = (l - ap) qi{x,t), 


( 110 ) 


where the hat denotes the corresponding solutions for a = 0 case, i.e., for Brownian particles. In the 
rest of this paper we shall follow the same notation. 

To derive the second cumulant we need to determine Yi, Ai and i? 2 - Combining the series expansion 
for Y with the definition (51) results in 


1 

Yi = - dx{qi{x,T) - qi{x,0)). 
P Jo 


Thanks to (110) we have 


Yi = (l-ap) Yi. 

To compute we start with the general formula ( [2^ which in the present case becomes 

pT 


( 111 ) 

( 112 ) 


R = 


pq(x,0) 

dx / dr 


1 


r(l — ar) 


pi pOO 

{q{x,0)-r)+ dt dx q{l - aq) {dj;pY. (113) 

Jo J — oo 


At second order we get 


i?2 = 


2p(l- 


1 f'CCi i>T nOO 

-r/ dx {qi{x,0)y + p{l - ap) dt dx {d^piY (114) 

^P) J— OO Jo J— oo 


which in conjunction with (110) yields 


i?2 = (l - ap)i?2- 


The above relations (112) and (115) also establish 

Ai = Ai. 


Using (112) and (116) we transform ( |102[ ) into 

(A^= (1 - ap) (aI 


The variance for the Brownian case is derived in (65) which then leads to 


= 2 

p a/tt 


(115) 


(116) 


(117) 


(118) 


This expression agrees with the general formula for the variance derived earlier in Section [4.2[ validating 
the approach we used here. 

The variance can also be determined by directly solving for qi and pi. An extension of this approach 
to the fourth cumulant requires explicit solution of the optimal fields up to the third order which is 
a very tedious task. The alternative approach based on the mapping to a = 0 case, as demonstrated 
above for the variance, considerably simplifies computations, so we adopt this approach in the following 
derivation of the fourth cumulant. Before proceeding, we write explicit formulas Ai and Yi. Using (56) 


which describes the case of a = 0, we compute Ai = p leading [due to (116)] to 

Ai = p. 


(119) 


Fi = (1 - ap) Vt. 


Similarly we derive 


(120) 
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The fourth cumulant 

In the second order in B we have 

{dt + dj:oo)p 2 = -(l - 2 ap) (do^pi)^, 

{dt - da:x)q 2 = -2/3(1 “ ap)“ 2(l - 2 ap)d^{qida:Pi), 

whereas the corresponding equations for the third order are 

{dt + dxa:)p:i = -2(1 - 2 ap){dxPi){dxP 2 ) + 2 aqi{dxPiY, 


{dt dxx^qs — dx 


2 p{l - ap)dxP3 + 2(1 - 2 ap)qidxP 2 + 2((l - 2 ap)q 2 - a{qi)'^)dxPi 


These equation are derived from the Hamilton equations (18|-(19). The corresponding boundary con¬ 
ditions follow from (104)-(105): 


P2{x,T) = -Yi S{x), 

q 2 {x, 0 ) = p{l - ap)p 2 {x, 0 ) + — — [pi{x,T) -pi(a;,0)]^ , 

P 3 ix,T) = d'{x), 

q 3 {x, 0 ) = p{l - ap)p 3 {x, 0 ) + \ gi (a;, 0 ) 92 ( 2 :, 0 ) - ^ ^ [ 9 i(x, 0 )]^. 


P(1 - ap) 


3[p(l - ap)f 


From the above formulas, one can verify that the solutions at the second order are related to the 
corresponding solutions for a = 0 case by a simple transformation: 


P 2 {x,t) = (1 - 2ap) P 2 {x,t) 
q 2 {x,t) = (1 - ap){l - 2 ap) 92 ( 2 :,!) 


_ ’ 


<?, _■ 

^i^(l-2ap) 


( 121 ) 

( 122 ) 


On the right-hand side we write Fi —>■ implying that Yi should be replaced by in the 

corresponding solution for a = 0 case obtained by treating Yi as parameter. 

A similar relation can be derived for the third order terms: 


+ a{l — ap)u(x,t), 

(l-2ap) 






1“:' (l-2ap) 


P3ix,t) = (1 - 2 ap) P3{x,t) 

93 ( 2 ;,!) = (1 - ap){l - 2 apY %{x,t) 

Here u{x, t) is the solution of 

{dt + dxx)u = 2 9 i {dxPi) 

subject to u{x, T) = 0, while h{x, t) is the solution of 

{dt - dxx)h =-2 p dxxU Y2dx [ 9 ? dxPi] 


+ a{\ — apYh{x^t). 


(123) 

(124) 

(125) 

(126) 


1 ^ 3 

subject to h{x,0) = p u{x,0) — — [ 91 ( 0 ;, 0)] . Note that both u{x,t) and h{x,t) do not depend on a. 

These functions also appear in the analysis of the statistics of time integrated current in the symmetric 
exclusion process on an infinite line (see Appendix [D| . 

In order to compute the fourth cumulant we must find i ?4 and Y 3 . Combining the series expansion 


(106|-(107) and (113) we get 


/ oo nT noo 

dxUi{x)+ / dt dxVi[x,t) 

-00 0 J — 00 
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with 


C/4 = 


[q2{x, 0)]2 + 2^1 (x, 0)q3{x, 0) 


+ 


2p(l - ap) 

(1 - 2 ap)^ 


(1 2 ap) .^q,^(2.^o))^g2(a;,0) 


2p2(l — ap)2 


(91(3;, 0)) 


(gi(x,0))" 


12p3(l - ap)3 12p2(l_ap)2 

V 4 =p{l-ap) {dxP 2 Y + 2{da;Pi) (d^Ps) - aqfidxPiY 
+ (1 - 2 ap) q 2 {dxPiy + 2 qi{dxPi) {dxP 2 ) 

We now again express R 4 through the corresponding solution for a = 0 case: 


i?4 = (1 — ap) (1 — 2 ap)‘ 


Ra 


Ya- 


Xi 

l-2a 


+ p^a{l — apY 21. 


(127) 


The last term is given by 


1 = 


dx 


dt 


qi{x,Q)\ (h{x,Q)\ ^ I (qi{x,Q) 
"12 


P 


dx 


‘ 2 .{dxPi) 


dxU 

P 


{dxPi) 


This only involves the solution for the a = 0 case. 

The second quantity required to compute the fourth cumulant is Y 3 . Combining the series expansion 
(991 and (511, and using the relation of the optimal fields to their counterparts for a = 0, we can express 

h{x,T) h{x,0) 


Y3 in terms of the hat variables: 


nc 

¥3 = [l-ap){\-2ap)^ % ^ +ap{l-apY 

^1-*- l-2op Jo 


dx 


(128) 


Combining this equation with (1271, we find that Y 3 — R 4 /p which appears in the fourth cumulant 
(1031 can be written as 


- — + Vs = (1 - ap){l - 2 apf' 
P 


-^+Y3 

P 


+ ap{l — apY 


-1 + 


dx I 


h{x,T) h{x,0) 


P" 


(129) 


The term inside the square brackets in the bottom line of (|129|) simplifies to 

h{x,T) ft,(x,0)^ 


-I- 


dx 


P" 


1 /8 - 6v^ 

4! 4 /^ 


Vf. 


(130) 


This identity was verified numerically using Mathematica. We also obtained an analytical proof by 
comparing to the analysis of the current in the SEP with uniform initial profile (see Appendix [D| . 

On the other hand, from the exact solution (54) for the a = 0 case we get (by substituting Y{B) 
from (99) in (54) and expanding in powers of B) 


-^+Y - ^ 

p ^ i2v^ 6 ' dVrv 




(131) 
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This allows us to simplify the term inside the square brackets in the top line of (129) 


-^ + Y, 

p 


Vi 




120F 6(1 - 2ap) 4VtV(1 - 2ap)^' 


Substituting Yi from (120), we obtain 


1 

7t 


P 


(1 — ap) 1 (1 — apY 1 


-■ l-2a 


12y^ (1 — 2ap) 'Ay/tt (1 — 2apY 


After all this work Eq. ( 129|) becomes 

i?4 


1 


Y. 


= (1 - ap){l - 2apf - (1 - apf{^ - 2ap) 

IzyTT oyTT 

3 1 ^ otp{l — apY f — 6 -\/ 2 ^ 


+ + 


4! 


(132) 


(133) 


(134) 


Plugging (134) and Ai = p into (103) we obtain 


(Yl = 


^ 1 — ^4 — (8 — 3V2^ap^ (l — ap) H-(l — ap)‘ 


2Vt 


(135) 


cumulant for the SEP 


Specifying (|135|) to a = 0 we get back the result (|66|), whereas setting a = 1 we arrive at the fourth 

2VT 


(^tI = 


1-p 


1 - ^4 - (8 - AV2)p^ (l - p) + ^ (l - pY 


(136) 


This result, announced in Ref. [^, i s valid in the annealed settingj^ When p —1, this expression 
matches the result derived in Ref. |41|. 

Remark: The fourth cumulant in the quenched setting can be calculated along similar lines to yield 


2(1 — 2apY arctan — 1^ + Q!p(l — ap) (4 — aV2) Vt. 




For Brownian particles, a = 0, we recover Eq. (|64|. Setting a = 1, we obtain the fourth cumulant for 
the SEP in the quenched case: 


(y )f = (1 - p) 


(l-2p)V9 


— arctan 

TT 


( 2^2 


-11 + (4 - 3^2) 


-Vf. 

TT 


6 A microscopic derivation for the Brownian point particles. 


The microscopic problem was first studied by Harris 17 , who derived an exact formula for the variance 


of the tagged particle. The analysis used the fact that trajectories of the particles are related to the 
trajectories of non-interacting particles with an exchange of particle index to keep the ordering same 
(see 77 for a recent reference). 

There is an equivalent description of the problem in terms of the phase space trajectories 20p5p0 . 
Consider 2n + 1 point particles diffusing on a one-dimensional line. The only interaction between 
particles is the hard-core repulsion which preserves the order of the particles. The particles are indexed 
by {—n, —n + 1, • • • , n}. The central particle is set to be the tagged particle. Let Y = {y-n, ■ ■ ■ , Vn}, 
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Fig. 6 A schematic representation of a sample trajectory of two one-dimensional Brownian point particles on 
the coordinate plane (xi,X 2 )- Due to the hard-core repulsion between particles the motion is confined in the 
domain X 2 > xi, whose boundary is denoted by the diagonal line. 


with —L < j/_„ < • • • < j/„ < L, be the positions at time t = 0. The positions at time t are denoted by 

The evolution of the particles can be described by diffusion in a 2n -I- 1 dimensional space confined 
to the chamber x_„ < X-n+i < • • • < x„ (see Figurefor a schematic). The probability Pt(X|Y) of 
the particle position follows a diffusion equation 


n 

dtPt{X\Y) = ^ dlPtiX\Y) 

i——n 


(137) 


where the diffusion coefficient is again set to unity. The single-file constraint is implemented by a 
reflecting boundary condition along the boundary of the chamber: 


dxiPt{X\Y) = 9a;.^jPt(X|Y) at Xi = for all — n < i < n. 


The solution to Eqs. (137)-(138) can be written as 





dkny^^ ^ 

a j=-n 


(138) 


(139) 


where cr is the permutation operator acting on the indexes of the particles. The above solution is 
applicable only within the chamber x_„ < x_„+i < • • • < x„. Equation (139) is one of the simplest 
examples of the Bethe ansatz. 

The Gaussian integrals in (1391 can be evaluated yielding 


1 2 

Pt(X|Y) = X! n gt{x\y) = . (140) 

CT j=-n ^ 


‘ We also derived (1361 by a straightforward perturbative calculation, without introducing the interpolating 
parameter a. 















25 


6.1 Probability of the tagged particle position 

The central particle {j = 0) is chosen to be the tagged particle. Without loss of generality, we assume 
that the tagged particle starts at the origin, yo = 0. The probability of finding the tagged particle at 
xo = X at time t is 


Probt(a;|Y) 



dX2 


dXn -Pt(X|Y). 


This function has been studied in great details by Rodenbeck et al. for the annealed initial condi¬ 
tion. Their results allow one to extract the large deviation function. Here we present a short alternative 
derivation, and then we determine the large deviation in the quenched case. 

Let 


G{x,yi,--- ,yn) 



(141) 


In terms of this function 


Probt(x|Y) = ^gt(x|y^(o)) G (-X,• • • ,-ya(-n)) G (x, • ,2/a(n)) ■ 

(7 


Let us rewrite the above formula as 

n 

Probt(x|Y) = ^ gt{x\yk)Ak[x,Yk], (142) 

k——n 

where we grouped the terms with same values of cr(0) = k. Further, Y^ denotes the subset of Y 
excluding the A:-th element y^. The amplitude is 


Ak [x, Yfc] — G ( X, 2 /(T),(- 2 )! ■ ■ ■ I yak{-n)) G {x,ycr^.{l),yak{‘2.)i ‘ ‘ ‘ TUakin)) j (143) 

<Tk 


where the permutation operator Uk denotes permutations acting on the subset Y^. 

Equation (143) can be simplified using an identity (174), proved in the Appendix]^ resulting in 


Afc [x, Yfc] = ^ ^ ^ ^ n 


E-n Cfc-i Ek+l I 

where eg are binary variables taking values ±1. 


e^k 


Q- 


X-ye 

y/U 


(144) 


6.2 Comparison with a random field Ising model 


The appearance of binary variables Cj in Eq. (144) suggests to seek a connection to an Ising model 
with non-interacting spins Cj. Such a connection indeed exists. To see it we use an identity 


V a/R ) cosh(/ij) 


where hj is the magnetic held acting on the spin e^, dehned as 


hj = arctanh 


erf 


Vj 


V4t 


This leads to an expression for the amplitude 

Ak [x, Yfc] = 


Z\\ik 


nj^fc2cosh(/ij)’ 


(145) 


(146) 


(147) 
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where Z [hfc] is the partition function of an Ising chain of 2n spins in the ensemble with total magne¬ 
tization zero, with = {h-n, ■ ■ ■, hk+i ,..., hn). The partition function is defined as 




[^'=1 = Z! ■ ■ ■ Z! Z! ■ ■ ■ Z! 


e „,0 


€fc-lCfe+l 






(148) 


The spins do not interact with each other. In terms of the partition function (148) the probability of 
tagged particle position can be written as 


Prob((a;|Y) 



9t{x\yk) 

nj^fc2cosh(/ij) 


Z [hfc]. 


(149) 


In this formulation in terms of the Ising model, the quenched case where the initial particle positions 
Uj are fixed corresponds to a quenched magnetic field. The annealed case where the yj are fluctuating 
corresponds to the fluctuating magnetic field. It is well known that the properties of a random field 
Ising model is different in the quenched and in the annealed ensemble. Then, the difference in the 
statistics of the tagged particle position between the quenched and the annealed initial condition can 
be related to the non-eqnivalence of the random field Ising model in the two ensembles. 


6.3 Normal diffusion 


There is a finite number of particles, 2n -I- 1 in our setting, diffusing on an infinite line. Therefore the 
tagged particle has a normal diffusion in the large time limit. This large time is set by t ^ since 
initially particles are distributed within an interval [—L,L]. 


In this limit, the leading behavior of the probability given by (149) comes from replacing yj = 0 for 
all j. This corresponds to the case where the magnetic field hj = h = arctanh[erf (a;/-\/4t)] is uniform. 

The corresponding partition function is Z[hfc] = Note that it is independent of h because the 

total magnetization vanishes. 


Substituting yj = 0 for all j into (149) we obtain 


Probt(a:|Y) ~ gt(x|0) 




{2n + l)l 


7 ,! r).\ 


(150) 


The result does not depend on the initial positions Y of the particles, that is whether the initial state 
is annealed or quenched. For large n the leading dependence on n becomes 


Probt (x = ^V^') — —^ (151) 

with '(/i(^) = ln[erfc(^)erfc(—^)]. At large times, the distribution is expected to be Gaussian. This can 
be confirmed by expanding '0(0 small ^ to give 

with the self-diffusion constant T> = ^. The self-diffusion constant decreases as n increases indicating 
the sub-diffusive behavior observed in the macroscopic calculation. 


6.4 Sub-diffusion 

This limit is defined by the number of particles n —>■ oo and L —^ oo keeping the density n/L = p 
constant. The central (tagged) particle is caged by infinitely many particles. From the hydrodynamic 
result in earlier sections it is expected that in this limit the probability of the tagged particle position 
has the large deviation form Probt(a;) x exp [—-s/it (j){x/y/^)\. We now derive 0(0 for quenched and 
annealed settings, and compare with the results from the hydrodynamic approach. 
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Quenched Case 


As a quenched initial state we choose the equidistant one: All particles at t = 0 are placed determin¬ 
istically with separation p~^ between adjacent particles. Then the initial position of the particle is 
ye = p~^i. We write the probability given by Eq. (142) as 


with 


Prob((a;^ 

) = Po{x) + P-{x) + P+{x), 

(153) 

Poix) 

= gt{x\0)Ao [x, Yo], 

(154) 

P-{x) 

—n 

= X] 9 t{x\yk)Ak [x,Yk ], 

(155) 


k=-l 


P+{x) 

= '^9t{x\yk)Ak [x,Yk]. 

(156) 


k=l 

All the three terms Pqj P± have the same asymptotic large deviation form at large t limit. They 
only differ in the sub-leading terms in t. The advantage of using the representation (153) is that it is 
much simpler to analyze the term Pq and to extract the large deviation function. We derive the large 
deviation function using a saddle point analysis of Pq. 

We use the expression for Aq from (144) which leads to 


Po{x) = gt{x\0)J2---J2J2---J2 


e-n e-1 ei tn t^O 


ie- 


X — p 


-It 




(157) 


We have used the quenched initial particle position ye = p The Kronecker delta function can be 
replaced by an integral, 


Po{x) =9tix\0) j dB [' 






1 


-erfcl ee 


X — p 

Vif 


where B is the integration variable and the factor 1/2 is for later convenience. In this form, the binary 
variables ee are decoupled and their sums can be performed rather easily, leading to 

Po{x) =gt{x\Q) J dB exp ^ log |l-k (e^ - l) ^erfc^ ^^^ | 


Taking the limit n —> oo with uniform density of particles p and replacing the summation over i by an 
integral, yields. 


Petix) = 


gtix\0) j dB exp p J d2:log |l-f (e^ - l) ierfc^^^^^ | 

To compute the large deviation function we write ^ = xfyf^ and recast the above formula into 

Po(CV^) = ^ J dB exp p-\/4t J dzlog |l-I-(e^ — l) ^erfc(z)| 

+pV^ J dzlog |l-I-(e“'® - l) ^erfc( 2 ;)| . 
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At large t the integral over B is dominated by the saddle point leading to a parametric solution of the 
large deviation function 

= =-pI cizlog|l+(e^-l)ierfc(.)} 

—p J dzlog |l + (e“'® — l) ierfc(z)| (158) 

with B determined from the optimization requirement = 0. These results are identical to 

(47)-(48) obtained using the MFT. 

Annealed Case 

In this case, the initial positions of the particles are uniformly distributed in the interval [—L,L\. 
Averaging over the initial positions the probability of the tagged particle position can be expressed as 

,1 rO ry-i 


Prob((a:) = 




dy- 


2 . . . 


-L 


tr 


dy-r. 


7T / ■ ■ 

fL 

• / dyn 

P Jo Jyi 

dyn-l 


Probt(a:|Y), 


(159) 


where Probt(a;|Y) is given in (142). As in the quenched case, we write Probt(a:) as a sum of three 
terms. All terms have the same large deviation form. Hence it is sufficient to analyze 


Probt(a:) - ^ 0 ( 2 ;) = gt{x\0) 


L” 


dy_ 


dy- 


2 • • • 


dy-r 


f-L 


'-L 



fL 

■ / dyn 

^ Jo jyi 

dyn-1 




e_i El e„ 

The symmetry of the integrand allows us to re-write it as 




Cf- 


X- ye 

V4t 


Po(2;) = <7(a;|0) ^ ^ ^ ^ 


n 7 / 9 


U=1 


X- ye 


The large deviation function can be derived from the above expression using a saddle point approx¬ 
imation. Replacing the delta function by an integral representation we get 


Poix) = 5t(2:|0) / dB 


HE' 

J=-l €t 


eeB/2 


1 

L 


’-L 


dye^eAc( eg 


HE 

e=l f-e. 




1 


dy^-erfc ee 


X- ye 

Vu 

X- ye 

y/M 


Completing the summation over the binary variables ee we get 


Po(2:) = 5t(a;|0) / dB 


1 


l + {e^ -l) - / dy-erfc 


l-L 


x-y 

y/At 
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where we have replaced yi by y, as the former reduces to a dummy variable after the summation. 
In the limit of n —>■ oo and L —> oo keeping njL = p fixed, the above expression results in 


P^{x) = 


gt{x\Q) J dBexp p (e® - l) J + p {e ^ - l) dyierfc^^^^ . 

In terms of rescaled coordinates ^ = xj^/ii and z = y/V^ the above can be written as 

Po(g-\/4^) = ^ J dBexp —l) J dz erfc(^ — z) + (e“^ — l) J dzerfc(z —^)| . 

At large t, the integral is dominated by the saddle point. This leads to the large deviation function 

(e---l) 


^ ,. logPo(Cv^) (e^-1) 

MO - - lim -A=-^ = -p- -z—^ 

t->oo 2 




with B determined from the saddle point condition 

d(l)A{0 

dB 


dz erfc(z) — p- 


= 0 . 


dz erfc(z), (160) 




(161) 


To show the equivalence with the MFT result (59) we use (1611 to deduce 


dz erfc(z) 


6 ^^ = 




dz erfc(z) 


(162) 


Combining this with (160) one recovers the large deviation function (59). 


7 Summary 


We studied the full statistics of the displacement of the tagged particle in single-file diffusion. Our 
analysis mostly relies on the macroscopic fluctuation theory (MFT). We found the full solution for the 
simplest single-file system composed of impenetrable Brownian particles (we computed the optimal 
paths, the large deviation function, etc.). This single-file system is also amenable to an exact microscopic 
analysis. In the annealed case, the large deviation function was originally computed in 20 and we 


presented another (shorter) derivation; the quenched case has not been studied before. The predictions 
based on the exact analyses and those which were derived using the MFT match for the single-file 
system composed of Brownian particles. 

Single-file systems with arbitrary transport coefficients cannot be solved exactly. Yet such general 
systems are tractable perturbatively. Specifically, we performed an expansion in the powers of the 
fugacity and derived an exact formula for the variance of position of the tagged particle valid in the 
general case of arbitrary D{p) and (j{p). The forth cumulant can also be computed for a fairly general 
class of single-file systems, e.g., for systems with constant diffusion coefficient and quadratic (in density) 
mobility. The well-known example of such system is the SEP, and we computed the forth cumulant of 
the tagged particle position for the SEP. 

Our study can be extended in many directions. One should be able to use the MFT to calculate 
statistical properties of the tagged particle trajectory (the two-time correlation functions 69 , the first 


passage times etc.). It would also be interesting to investigate the influence of a bias in the simple 
exclusion process. If one consid ers an asymmetrically hopping tagged particle in a bath of SEP, all the 
cumulants scale as \/T 
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The full statistics has been recently determined in the high density 
limit 41 when the leading asymptotic can be extracted by treating the vacancies as independent 


random walkers; for the finite density even the variance remains unknown. If all particles are biased, 
even the scale of fluctuations depend o n th e setting: In the annealed case, the position fluctuation 

; in the quenched case, the exponent changes from 1/2 
lether a hydrodynamic approach based on the MFT can 


13 


of the tagged particle grows as OT [9 
to 1/3 [^ [11|[T4{[83}]8^ . It is not known w 
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capture these behaviors and if the Tracy-Widom distribution, which is ubiquitous in these problems 
(see e.g. [^), can be retrieved as a solution of some optimal path equations. 

From a more general point of view, it would be interesting to give a physical interpretation of the 
conjugate field that appears in the MFT optimal path equations and to classify the one-dimensional 
gases that could lead to classically integrable partial differential equations, that can be solved using 
the inverse scattering method 
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A Functional derivative of Y\q] 


The tagged particle position is a functional Y[q\ of the initial q{x, 0) and the final q(x, T) density profile, defined 
by the relation 


|•Y[q\ na 

/ dx q{x,T) = / 
Jo Jo 


dx 


q{x,T) - q{x,0) 


(163) 


which comes from equation This way of writing makes both the integra ls c onvergent. 

The functional derivatives appearing in the boundary conditions ( |15[ |- |l7| can be easily computed from 
the definition above. To proceed, consider a small variation in the final density q{x, T) —>■ q{x, T) + Sq{x, T) 
leading to a change Y ^ Y + 5Y. Gorresponding to this variation the above equation becomes 


L 


Y+5Y 


dx 


q{x,T) + Sq{x,T) 


f 


dx 


q{x,T) - q{x,0) 


roo 

- / dx 6q{x^ T). 
Jo 


Using the formula ( |163[ l this reduces to, 

rY+5Y 


p Y +0 Y POO 

/ dx q{x^T) = / dx 6q{x^T). 
Jy Jy+5y 


(164) 


Assuming that the change 5Y is small for a small variation Sq{x, T) and keeping only the linear terms, the 
above leads to 


SY 



'0{x-Yy 
. q{Y,T) _ 


5q{x,T). 


(165) 


Hence, we have 


SY _ 0{x-Y) 
Sq{x,T) q{Y,T) 


The functional derivative with respect to the initial density q{x, 0) is similarly derived: 


(166) 


SY _ 0(x) 

Sq(x,0) q{Y,T)' 


(167) 


B Inequality between cumulants in the annealed and quenched cases 


Because of the fluctuations in the initial state, the tagged particle is expected to have larger displacements in 
the annealed case than in the quenched case. More generally for any single file system 


/i^(A) > ms(A) 


(168) 


which implies ^ 

To derive (|168l) we recall that the two initial settings differ by how the average over initial state is taken: 


/i^(A) = log<;^e^^^^ 

MS (A) = ^logll^e^-^'^^ 


evolution+initial 
evolution^ initial 


(169) 

(170) 


Here the subscript evolution denotes average over stochastic evolution of the system and t he s ubscript initial 


denotes average over initial state. E quati on (1681 then follows from the Jensen inequality 87 because log is 


a concave function. The inequality (1681 implies the opposite relation for the large deviation functions and 


therefore explains why the large deviation function in the quenched case exceeds the large deviation function 
in the annealed case (Fig. |^. 
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C Mobility cr{p) 


We first derive o'(p) for the Brownian particles with hard core repulsion. Let there are n number of particles 
at equilibrium within an interval of length L. All particle positions within the interval are equally probable. 
This leads to the canonical free energy density 



In the limit n —>■ oo and L —>■ oo with finite n/L — p, the above formula becomes /(p) = plnp — p. Using it 
together with the fluctuation-dissipation relation one gets a{p) = = 2 p. 

Consider now the SEP with n particles on the ring of L sites. In the equilibrium, all configurations are 
equally probable which leads to the free energy density 



Using Stirling’s approximation we get /(p) = plnp+(l — p) ln(l — p). Using the fluctuation-dissipation relation 
we recover the well-known formula (j(p) = 2p(l — p) for the mobility in the case of SEP. 


D Derivation of the integral (130) 


The result in ( |130[ ) can be proved by cornparing with the analysis of the time integrated current in a symmetric 
exclusion process on an infinite line [M 73 . Consider an annealed initial state at uniform density p. Let Qt 
be the time integrated current in a time window [0, T] through the site at origin. The current is related to the 
hydrodynamic density profile by 


Qt = 


/■ 

^0 


dx 


p{x,T) - p{x,G) 


On an average the current is zero because of the uniform initial density. However, all the even cumulants 
are non-zero. The analysis for the cumulant generating function of Qt can be easily formulated in terms of the 
macroscopic fluctuation theory [59| . To analyze the variational formulation we take a series expansion method, 
similar to the one presented in Section^ This way the fourth cumulant of current for the symme tric e xclusion 
process can be related to that for non-mteracting particles. The relation is similar to the one in (|129||: 


Qt) = (1 — ap)(l — 2 ap)^ (Qr) -|-24 qp^( 1 — ap)^ 


-1 + 


r 


dx 


h{x,T) h{x,0) 


(171) 


where I and h are defined in section We followed the same convention of denoting the non-interacting case 
by hat variables. 

Using the Bethe ansatz Derrida and Ger sche nfeld derived an exact expression for all the cumulants of the 
current for the SEP in the annealed setting 73 . The fourth cumulant in the annealed case can be extracted 
from their work: 

l-3%/2p(l-p) 


— p(l — p) 


(172) 


On the other hand t he c umulant for non-interacting particles corresponding to a = 0 can be derived by a direct 
counting argument leading to 


Qt') — P —j=- 

/ c -/tT 


(173) 


Making use of these results we establish (130l. 


E Function G 


Here we prove an identity for the function G defined in equation (1411. Let 7 be a permutation of the set 
17 = {y-n, ■ ■ ■ , y-i,yi, ■ ■ ■ , yn} of 2n elements. Note that the zeroth element has been excluded for convenience. 
The function G has the property that 


^ ^ ( ^5 y ^{ — n ) i ''' ) 2/7( —2}) 2/7( — 1}) (^1 1/7(1) 5 1/7(2) 5 ' ■ ' ^ y ^( n )) 

7 


e^o 


X-ye 


(174) 
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The identity holds for any set of 2n elements, with n being a positive integer, and the binary variables ej = ±1 
for any j. 

First, we want to prove 


where u is a permutation of the set {yi, ■ • • , i/„}. Using the definition of G in ( |141[ ) we get 

POO POO POO 

= ^ / dzi dz2 ■ ■ ■ dz^Yl gt {zj\y^(j)) . 

U u j = l 


(175) 


As Zj^s are dummy variables, the expression can be re-written as 


^ ^ yiy(i)' ’' 1 2/i^(n)) — 


POO POO P 

I] / dz.(i) / da,(2) ■■■ 


dz, 


v(n) 




i=i 


The range of integration denoted inside the square bracket is equal to the integration over the entire volume 
of all Zj > X, leading to 


^ ^ G(x, ■,yu{n)) 


n / 


dzj gt izj\yj 


n 

=n 




X - yj 

y/u 


(176) 


proving (1751. 

In the next step we prove the identity (1741 using (1761. In (1741 the summation is over all permutations 
of the set Q. Any ordered set which is generated by a permutation of O can also be generated by suitably 
constructing two subsets of n elements each, and subsequently taking permutations within the elements of the 
subsets. To elaborate, let us divide the ordered set Q = {y-n, ■ ■ ■ ,y-i,yi, ■ ■ ■ ,yn} into two ordered subsets 
ljJl = {y-n, ■ ■ ■ , y-i} and ujr = {y\, • • • , yn}- The number of ordered sets generated by permutations of Q is 
(2n)!, all of which can also be generated by combination of the following steps. 

1. The simplest are permutations strictly within ujl and ujr. There are n\n\ such permutations. 

2. One can construct another pair of ordered subsets uj'r and uj'r by exchanging any two elements between 
UJL and ujR. There are pair of subsets that can be generated this way. Overall this produces n!n!(")^ 
permutations. 

3. We can further exchange two elements from ujr with two elements from ujr, and then consider permutations 

within those subsets. This generates permutations. 

4. Continuing this process we can generate all permutations O. 

As a useful check, one can compute all the permutations generated by the above procedure and find that it is 
equal to the total number of permutations: 




1 + 


n \ In 

1 + 2 


+ ■ 


= (2n)! 


The above decomposition of permutations simplifies the summation on the left hand side in (174|. For 
example, consider the ordered sets generated in step 1 of the decomposition where the ordered subsets are 
UioJl and V 2 <^r, with Vi and U 2 being permutation operators. Summing over all permutations we have 


G( X, yu^(^ — n),''' , yiyii — l)) ^ ) G{x, l/iy 2 (l); q ( 


i=i 


X - y-j \ 1 


-erfc 


x-yj 


where we used the relation ( |175| ). We now introduce binary variables {e-„, ■ • ■ , e_i, ei, • ■ • , e„}, each of which 
take values ±1. Then the above formula can be re-written as 

y G(—a;, — , ■ • ■ , — (_!)) y ( G{x, 2 / 1 / 2 (i ),''' , y^ 2 (n)) — 9 7zT~ ) ’ 

I/, I/O ^ V V4r / 


with a configuration of the binary variables: £j = — 1 for j < 0 and ej = 1 for j > 0. 

The summation over configurations generated in the step 2 is performed in a similar way. In this, the 
subsets are generated by exchange of one element each from cur and ujr. For example, consider the (—n)th 
element from ujr and 1st element from ujr have been exchanged. The summation over configurations generated 
by permutation within these subsets yields a simple expression 

y G(—X, —2/1/1 (1) , — 2/1/1 (-Ti-l-l) ■ ■ ■ 5 —2/1/1 (-1)) y G{x, y,^2{ — n), yv2 ( 2 )) ■ ■ ■ I 2/1/2(11)) — I Ej / , ) , 

1/1 1/2 j = -n,j^0^ ^ V 4 il£/ 

with e_„ = 1, El = — 1 and other binary variables are the same as before. Fror n the se examples the pattern 
emerges, namely the above decomposition of permutations leads to the identity (174l. 





















33 


References 

1. Chou T and Lohse D 1999 Phys. Rev. Lett. 82 3552 

2. Karger J and Ruthven D 1992 Diffusion in zeolites and other microporous solids (New York: Wiley) 

3. Das A, Jayanthi S, Deepak H S M V, Ramanathan K V, Kumar A, Dasgupta C and Sood A K 2010 ACS 
Nano 4 1687 

4. Gupta S, Rosso A and Texier C 2013 Phys. Rev. Lett. Ill 210601 

5. Li G W, Berg O G and Elf J 2009 Nature Physics 5 294 

6 . Hodgkin A L and Keynes R D 1955 Journal of physiology 128 61 

7. Richards P M 1977 Phys. Rev. B 16 1393 

8 . De Masi A and Ferrari P A 1985 J. Stat. Phys. 38 603 

9. Kipnis C 1986 Ann. Probab. 14 397 

10. Majumdar S N and Barma M 1991 Phys. Rev. B 44 5306 

11. van Beijeren H 1991 J. Stat. Phy. 63 47 

12. Ferrari P A and Fontes L R G 1996 J. Appl. Probab. 33 411 

13. Liggett T 2004 Interacting Particle Systems (New York: Springer) 

14. Gupta S, Majumdar S N, Godreche C and Barma M 2007 Phys. Rev. E 76 021112 

15. Percus J 1974 Phys. Rev. A 9 557 

16. Roy A, Narayan O, Dhar A and Sabhapandit S 2013 J. Stat. Phys. 150 851 

17. Harris T E 1965 J. Appl. Probab. 2 323 

18. Arratia R 1983 Ann. Probab. 11 362 

19. Sethuraman S and Varadhan S R S 2013 Annals of Probability 41 1461 

20. Rodenbeck C, Karger J and Hahn K 1998 Phys. Rev. E 57 4382 

21. Kollmann M 2003 Phys. Rev. Lett. 90 180602 

22. Imamura T and Sasamoto T 2007 J. Stat. Phys. 128 799 

23. Euan-Diaz E C, Misko V R, Peelers F M, Herrera-Velarde S and Castaneda Priego R 2012 Phys. Rev. E 
86 031123 

24. Flomenbom O and Taloni A 2008 EPL 83 20004 

25. Lizana L and Ambjornsson T 2008 Phys. Rev. Lett. 100 200601 

26. Ben-Naim E and Krapivsky P L 2009 Phys. Rev. Lett. 102 190602 

27. Manzi S J, Torrez Herrera J J and Pereyra V D 2012 Phys. Rev. E 86 021129 

28. Lizana L, Lomholt M A and Ambjornsson T 2014 Physica A 395 148 

29. Krapivsky P L, Mallick K and Sadhu T 2014 Phys. Rev. Lett. 113 078101 

30. Barkai E and Silbey R 2009 Phys. Rev. Lett. 102 050602 

31. Kukla V, Kornatowski J, Demuth D, Girnus I, Pfeifer H, Rees L V C, Schunk S, Unger K K and Karger 
J 1996 Science 272 702 

32. Wei Q H, Bechinger C and Leiderer P 2000 Science 287 625 

33. Lutz C, Kollmann M and Bechinger C 2004 Phys. Rev. Lett. 93 026001 

34. Lin B, Meron M, Cui B, Rice S A and Diamant H 2005 Phys. Rev. Lett. 94 216001 

35. Siems U, Kreuter C, Erbe A, Schwierz N, Sengupta S, Leiderer P and Nielaba P 2012 Scientific Reports 2 
1015 

36. Edwards S F and Wilkinson D R 1982 Proc. Royal Soc. A 381 17 

37. Kardar M, Parisi G and Zhang Y C 1986 Phys. Rev. Lett. 56 889 

38. Levitt D 1973 Phys. Rev. A 8 3050 

39. Alexander S and Pincus P 1978 Phys. Rev. B 18 2011 

40. Krapivsky P L, Redner S and Ben-Naim E 2010 A Kinetic View of Statistical Physics (Cambridge Uni¬ 
versity Press) 

41. Illien P, Benichou O, Mejia-Monasterio C, Oshanin G and Voituriez R 2013 Phys. Rev. Lett. Ill 038102 

42. Leibovich N and Barkai E 2013 Phys. Rev. E 88 032107 

43. Jara M and Landim C 2006 Annales de I’Institut Henri Poincare Probability and Statistics 42 567 

44. Derrida B 2007 J. Stat. Mech. P07023 

45. Touchette H 2009 Phys. Rep. 478 1 

46. Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2001 Phys. Rev. Lett. 87 040601 

47. Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2002 J. Stat. Phys. 107 635 

48. Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2005 Phys. Rev. Lett. 94 030601 

49. Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2009 J. Stat. Phys. 135 857 

50. Bertini L, De Sole A, Gabrielli D, Jona-Lasinio G and Landim C 2014 ArXiv e-prints 1404.6466 

51. Jordan A N, Sukhorukov E V and Pilgrim S 2004 J. Math. Phys. 45 4386 

52. Pilgram S, Jordan A N, Sukhorukov E V and Biittiker M 2003 Phys. Rev. Lett. 90 206801 

53. Derrida B 2011 J. Stat. Mech. P01030 

54. Bertini L, Sole A D, Gabrielli D, Jona-Lasinio G and Landim C 2006 J. Stat. Phys. 123 237 

55. Bodineau T, Derrida B and Lebowitz J 2010 J. Stat. Phys. 140 648 

56. Bodineau T, Derrida B, Lecomte V and van Wijland F 2008 J. Stat. Phys. 133 1013 

57. Tailleur J, Kurchan J and Lecomte V 2007 Phys. Rev. Lett. 99 150602 

58. Bunin G, Kafri Y and Podolsky D 2012 EPL 99 20002 

59. Derrida B and Gerschenfeld A 2009 J. Stat. Phys. 137 978 

60. Jona-Lasinio G 2010 Prog. Theor. Phys. 184 262 

61. Jona-Lasinio G 2014 J. Stat. Mech. P02004 

62. Krapivsky P L and Meerson B 2012 Phys. Rev. E 86 031106 

63. Meerson B and Sasorov P V 2013 J. Stat. Mech. P12011 



34 


64. Vilenkin A, Meerson B and Sasorov P V 2013 J. Stat. Mech. P06007 

65. Meerson B and Sasorov P V 2014 Phys. Rev. E 89 OIOIOI(R) 

66 . Meerson B, Vilenkin A and Krapivsky P L 2014 Phys. Rev. E 90 02210 

67. Hurtado P I, Espigares C P, del Pozo J J and Garrido P L 2014 J. Stat. Phys. 154 214 

68 . Krapivsky P L, Mallick K and Sadhu T 2015 J Phys A: Math Theor 48 015005 

69. Krapivsky P L, Mallick K and Sadhu T 2015 ArXiv 1505.01287 

70. Baek Y and Kafri Y 2015 ArXiv e-prints 1505.05796 

71. Spohn H 1991 Large Scale Dynamics of Interacting Particles (New York: Springer-Verlag) 

72. Rajesh R and Majumdar S N 2001 Phys. Rev. E 64 036103 

73. Derrida B and Gerschenfeld A 2009 J. Stat. Phys. 136 1 

74. Freidlin M and Wentzell A 1984 Random Perturbations of Dynamical Systems (New York: Springer-Verlag) 

75. Martin P, Siggia E and Rose H 1973 Phys. Rev. A 8 423 

76. De Dominicis C and Peliti L 1978 Phys. Rev. B 18 353 

77. Hegde G, Sabhapandit S and Dhar A 2014 Phys. Rev. Lett. 113 120601 

78. Elgart V and Kamenev A 2004 Phys. Rev. E 70 041106 

79. Prudnikov A P, Brychkov Yu A and Marichev O I 1986 Integrals and Series: Special Functions vol 2 (New 
York: Gordon & Breach Science Publishers) 

80. Kumar D 2008 Phys. Rev. E 78 021133 

81. Burlatsky S, Oshanin G, Mogutov A and Moreau M 1992 Phys. Lett. A 166 230 

82. Landim G, Olla S and Volchan S B 1998 Commun. Math. Phys. 192 287 

83. Gartner J and Presutti E 1990 Ann. Inst. Henri Poinc. 53 1 

84. Ferrari P A 1991 Ann. Inst. Henri Poinc. 55 637 

85. Alexander F J, Janowsky S A, Leibowitz J L and van Beijeren H 1993 Phys. Rev. E 47 403 

86 . Babelon O, Bernard D and Talon M 2003 Introduction to Classical Integrable Systems (Cambridge: Cam¬ 
bridge University Press) 

87. Gradshteyn I S and Ryzhik I M 2007 Table of Integrals, Series, and Products (London: Academic Press) 



